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Cosmic inflation provides a mechanism for generating the early density perturbations that seeded the large-scale 
structures we see today. Primordial non-Gaussianity is among the most promising of few observational tests 
of physics at this epoch. At present non-Gaussianity is best constrained by the cosmic microwave background, 
but in the near term large-scale structure data may be competitive so long as the effects of primordial non- 
Gaussianity can be modeled through the non-linear process of structure formation. We discuss recent work 
modeling effects of a few types of primordial non-Gaussianity on the large-scale halo clustering and the halo 
mass function. More specifically, we compare analytic and N-body results for two variants of the curvaton 
model of inflation: (i) a "tjvx" scenario in which the curvaton and inflaton contribute equally to the primordial 
curvature perturbation and (ii) a "gjvL" model where the usual quadratic /jvl term in the potential cancels, 
but a large cubic term remains. 



1. Introduction 

From the anisotropies in the cosmic microwave background (CMB) to the large-scale distribution of dark 
matter halos hosting galaxies, the universe appears rich with structure. A central goal of observational cosmology 
is to understand the period of inflation that is believed to have generated this structure [H-Q- Current CMB 
data confirms inflationary predictions for a spatially flat universe with primordial curvature perturbations drawn 
from a nearly scale-invariant power spectrum |5j. Nevertheless, distinguishing between microphysical models 
remains a challenge. 

Evidence for non-Gaussianity in the primordial perturbations could rule out large classes of inflationary 
models and shed light on the mechanism that generated the initial structure (see e.g. 6 L ZQ. At present the 
most stringent bounds on primordial non-Gaussianity come from CMB experiments [5|, |g, |9|, but data from 
galaxy surveys is increasingly competitive (see, e.g. [Iol - fl2| ). In these proceedings we give an overview of 
analytic and N-body results for the abundance and clustering of dark matter halos with local non-Gaussian 
initial conditions described by the parameters /nl, 3nl and tjvl- In fj2]we introduce three examples of non- 
Gaussian initial conditions. In S|3]we review analytic models of the impact of local non-Gaussian initial conditions 
on the abundance and clustering of dark matter halos. Comparison of the analytic models and results from 
N-body simulations is given 2) In Sj5]we summarize our results. The reader is referred to the original references 



iN-Doay 



for a more detailed discussion. 



2. Examples of local non-Gaussian statistics: J NL , g NL , r NL 

If the initial curvature 1 perturbations are homogeneous, isotropic and Gaussian their statistics are entirely 
characterized by the two point correlation function, or power spectrum 

($(k)$(lO) = (2vr) 3 ( 5 Dil . ac (k + k') J P*$(fc). (I) 

If $ has any non-zero odd Appoint function or even Appoint function that just isn't specified by the two-point 
function (i.e. ($(xi)$(x 2 ) . . . $(xjv)) ^ ($(xi)$(x 2 ))($(x 3 )$(x4)) . . . {§(x N -i)$(x N )} + permutations) then 
$ is non-Gausssian. 

While there are an abundance of inflationary scenarios producing a variety of non-Gaussian initial conditions 
(there are some organizing principles; see for example [13]) here we focus on initial curvature that can be written 
as a non-linear mapping of a Gaussian random field that is local in real space. These types of initial conditions 
can arise in the curvaton model (l8l.[l9|. For a general review of inflationary scenarios giving local non-Gaussian 
initial conditions see [201 - 



1 Following standard notation in studies of non-Gaussianity, we define <3? = where f is the gauge invariant primordial 
curvature. 
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2.1. /nl 

Consider defining the initial curvature as a Gaussian random field $ G , plus a small (0(/nl\/ ($ G ))) fractional 
perturbation [2l[ 

$Ar G (x) = <f G (x) + /nL ($ G (x) - ($ 2 G )) . (2) 

The new field <E>jv G (x) obeys non-Gaussian statistics. In particular it has a skewness (&%q) ~ G/nl^c) 2 , 
and kurtosis ($^ G ) C = ($ 4 NG ) - 3( c f )2 VG ) 2 — 48/f 2 JL (<I ) G ) 3 . The CMB constraints on this form of primordial 
non-Gaussianity are — 10 < /nl < 74 at 95% confidence |5|. 

Another non-Gaussian feature is the coupling between different physical scales. We can gain some insight into 
the mode coupling by splitting the Gaussian field into "short" and "long" wavelength pieces, <E> G = $ GjS + $ Gi /. 
The short wavelength fluctuations of the non-Gaussian field depend on both short and long wavelength modes 
of the Gaussian field, 

<f>NG.s = ®G,s + /nl ($ g , s - + 2/ NL $ G , s $ G ,i ■ (3) 

In particular, an observer on top of a long wavelength mode $ G ,; will see small-scale statistics that depend on 
the value of $ Gj z 

(^ngJI ~ (®ng, s ) (l + 4/NL<S> G ,* + 4/£ L $ G>i ) (4) 



<^G )S )|,~6/ NL ($ G) ^(l + 4/ NL $ G ,0 (5) 

where we have only kept terms up to 0(f^ L ). This coupling between short and long wavelength scales will be 
important to understand the impact of non-Gaussian initial conditions on the clustering of dark matter halos. 



2.2. # NL 

A variant on Eq. ^ is to consider a local mapping where the quadratic term vanishes, and the cubic term is 
important [H, Hl| 

$at G (x) = $ G (x) + <?nl (* G (x) - 3$ G (x)($ G }) . (6) 

With this mapping the skewness of $ng vanishes and the kurtosis is (&%q)c ~ 24<?nl('I ) g } 3 - The CMB limits 
-12.34 < c/nl/IO 5 < 15.58 at 95% confidence [|. 

In the (?nl model, the coupling of short and long scales in Eq. ^ gives a small-scale variance that depends 
on 

and a local skewness the varies linearly with &g,i 

<$WIz~ 18 5nl(<& g , s } 2 $ G ,z. (8) 

Comparing the leading terms in Eqs. ([5]) and (JSJ, we see that to an observer sitting on a long wavelength mode 
it appears that they live in a cosmology with /^{^(x) = SgNL'&G^lx)' 



2.3. r NL 

Another variation is to consider initial curvature which is the sum of two fields G and cr G which fluctuate 
independently ((0 G cr G ) = 0) but have proportional power spectra, P^/Paa = £ 2 = constant. Non-Gaussianity 
is generated by adding a term quadratic in a, 



$Ar G (x) = <Mx) + *G(X) + /NL(1 + t 2 ) 2 (4W ~ (4)) • (9) 
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The skewness in this model is given by ($% G ) ~ 6/nl( ( I > ^g) 2 , and the kurtosis is ($% G ) C ~ 48/ NL (l + 
£ 2 )(®ng) 3 = 48rAri($|, G ) 3 /(|) 2 where the factor of 6/5 is conventional. 2 Current CMB bounds on this 
parameter are —6000 < tjvl < 33,000 at 95% confidence 

In this tjvl model, there is a coupling between small-scale fluctuations in $>ng and the long- wavelength 
fluctuations in a G 

(^ s )L-(^ s )(l + 4/NL(l+C 2 )^+4/ 2 L (I+^ 2 ) 3 a 2 ; ,) . (10) 

In this two-field example tn j, > (§/nl) 2 - One may wonder whether this inequality is fundamentally related to 
the physics of inflation. [25| showed the inequality was true at tree- level using the SN formalism. In fact, this 
inequality can be interpreted as a positivity constraint that must be satisfied regardless of the mechanism that 
generated the perturbations. A formal proof is given in [l6[, but the argument can be understood as follows. 
From Eqs. (j4]) & (fT0| we can see that /nl is a measure of the large-scale correlation between the potential $ and 
the locally measured small-scale power (/nl ~ / (®f)(®1))- On the other hand, t^l is a measure of the 

large-scale variance in the small-scale power (tjvl ~ ($ 2 $ 2 ) c /($ 2 ) 2 ($ 2 )). The inequality tml > (|/nl) 2 then 
arises as the condition that the correlation coefficient between the small-scale power and $ must be between -1 
and 1. 



3. Impact of non-Gaussian initial conditions on large-scale structure: analytic 
predictions 

In the previous section we discussed non-Gaussianity in the initial curvature perturbation. To understand the 
impact of primordial non-Gaussianity on the abundance and clustering of dark matter halos, we need a model 
that relates halos to perturbations in the initial curvature. In linear theory, the matter density perturbation 6 
is simply related to the initial curvature, 

where T{k) is the transfer function and D(z) is the linear growth function. A simple prescription for halo 
formation is to smooth the linear density field on scale M 



(x) = I / d 3 x'<5(x - x') (12) 

v JV~M/p m 



and to model halos of mass > M as regions of the smoothed initial density field with <5m(x) > 5 C , where S c is the 
spherical collapse threshold. The probability distribution function (PDF) for fluctuations 8m can be written in 
terms of the cumulants for the smoothed density fluctuation 5m 



T M a M 



, M ^(<5 2 / ), K3 (M)^^, Ki (M) = ±M£, ... (13) 



where K3 oc /nl and K4 contains terms oc <7nl and cx tjvl- If o~m is non-Gaussian, an infinite number of 
cumulants are needed to completely specify the PDF. However, one can approximate the PDF with a finite set 
of cumulants and for our purposes a M , K3 and K4 are sufficient. 

In the next sections we use these elements to model the abundance and clustering of dark matter halos. In 
the plots throughout this paper we use the WMAP5+BAO+SN fiducial cosmology [26|]: baryon density Vl^ti 2 = 
0.0226, cold dark matter (CDM) density ft c h 2 = 0.114, Hubble parameter h = 0.70, spectral index n s = 0.961, 
optical depth r = 0.080, and power-law initial curvature power spectrum k 3 P^(k)/2iT 2 — A 2 (fc/fc p i v )" s_1 where 
A 2 = 2.42 x 10~ 9 and fc piv = 0.002 Mpc" 1 . 



2 Apparently because /nl is typically denned through Eq. J2} giving /jjl ~ ^{$ 3 } / {& 2 ) 2 , where $ = |f with ( the primordial 
curvature, but tjvl is typically defined in terms of £ as rjvL ~ (C 4 )c/(C 2 ) 2 1241 - 
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3.1. Halo abundance 

Dark matter halos form from rare positive fluctuations in the matter density. Press & Schechter [22} gave a 
simple analytic model for the abundance of dark matter halos with mass M in terms of the PDF for the linear 
matter density fluctuations smoothed on scale R = (3M/ 47rp ITl ) 1 /^ , 

n(M) = -2^JLp(5 M >5 c ,M) (14) 
where S c « 1.68 is the spherical collapse threshold and V(5m > S c , M) is the probability for Sm > S c . Using 

r°° i 

V G3bassian (6M>5 c ,M)= dS M -- e-'&Arj' (15) 

J Sc V 27R7 M 

in Eq. (fl4"|) gives the Press-Schechter mass function 1271 . which is known to disagree at the ~ 50% level with the 
mass function measured from N-body simulations [28|]. Nevertheless, the expression in Eq. (fT4"]) has long been 
used to predict the relative abundance of halos in a non-Gaussian cosmology to a Gaussian one [2{|. Typically, 
the full non-Gaussian PDF is not known, but approximate expressions for n(M) are obtained by truncating 
an asymptotic expansion (e.g. 30]) or Edgeworth series (e.g. [3lj - hereafter the "Edgeworth mass function") 
for the PDF and using Eq. (| 14[) . These expressions agree well with N-body simulations with /NL-type initial 
conditions provided the "modified" collapse threshold S' c w 1.42 is used (e.g. (32l|). 

In [l~3j], truncating the series for hi V(Sm > <5 C , M) rather than V(5m > <^c) was proposed, where the Edgeworth 
expression is used for the PDF. In the limit of small non-Gaussian corrections, the mass function obtained this 
way (which we call the "log-Edgeworth" mass function) agrees with the Edgeworth mass function, but the in 
high-mass limit where non-Gaussian corrections are important, the "log-Edgeworth" mass function is better 
behaved (see Figs. (JTJ) & ©). 



3.2. Scale-dependent halo bias from / NL and % L -type initial conditions 

In the previous section we discussed modeling dark matter halos as regions where the linear density field Sm 
exceeds 8 C . If a fluctuation 8m is sitting on top of a longer wavelength fluctuation in the density Si, then the 
local collapse threshold is adjusted to S c — Si, thus the halo abundance fluctuates with the density as 



Sn(x) dlnn(M,5 c — Si) 



dSi 



Si(x), (16) 

<5i=0 



where the 1 accounts for the fact that the Eulerian halo number density is increased by a factor of 1 + S with 
respect to the Lagrangian one given by n{M). 

As we've seen in non-Gaussian initial conditions can cause small scale statistics like the variance and 
skewness to be modulated by the long- wavelength potential For the /nl and <?nl initial conditions the small 
scale variance and skewness are modulated by Accounting for this modifies Eq. (fT6|) to 



5n , , Sinn, s d Inn , , 9 Inn , . , . 

— (x) = 1 + -jtt-Hx) + 4/nl 7 t 1 2"$j(x) + 3snl-t7— $j(x) ■ (17) 

n aoi oma M of nl 

In Fourier space, the density field is related to the early-time gravitational potential through a(k, z) which 
allows us to write 

Sn ( 4/ NL Sinn 3g NL 91n?i\ 

— (k) = [b+ — — + —7-J — r-oT— )^( k ) ( 18 ) 

n \ a{k,z)oma M a(k,z) c/jnl / 

= b fNM (k)Si(k) (19) 

where 6 = 1 + d\nn/dSi. We'll refer Eq. (fT5)l to as the "peak-background-split" (PBS) prediction for the 
scale-dependent non-Gaussian bias. 

For a mass function dependent only on the combination <5 c /(tm, rather than S c and o~m separately, we can 
write the a M derivative in terms of the constant bias b 

— (k) = (b + g^-i) + m . (20) 

n V a(k,z) a{k,z) 9/ NL J 
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The 5nl bias coefficient can be rewritten in terms of b, but the expression is more complicated [15| . The /nl 
dependent bias above was first written down in (33|, but the derivation given here follows that of [l(|. The 
expression for <?nl bias term is from 15], but see also (34l. l35j|. 

From Eq. (f2"U)) the large-scale halo-matter and halo-halo power spectra are given by 

Phm(k) = b fNht9Nh (k)P mm {k) and P hh (k) = b 2 fN ^ gNL (k)P mm (k) . (21) 

The analytic form for the /NL-dependence of -FW and Phh given by Eqs. ([20)) & (|2~T|) has been shown to be in 
excellent agreement with simulations (see e.g. 32, 33, 36]). Using the scale dependent halo bias, [l(| constrain 



-29 < /nl < 70 at 95% confidence. Using the fact that the scale-dependent bias from (?nl has the same 
fe-dependence (ex l/a(k, z)), [34j applied the /nl constraints from [l(| to limit —3.5 < <?nl/10 5 < 8.2. 

3.3. Stochasticity between halos and dark matter from r NL ^ (|/nl) 2 initial conditions 

For the two-field initial conditions given in ^2.31 the linear density field is determined by the sum of the 
potentials a + <j), through 

$(M)=a(M(0(k) + <7(k)) > (22) 
but the small scale power is modulated by 4/nl(1 + £, 2 )ci,g- So the number of halos fluctuates as 

on 

— (k) = b6i(k) + 2/ NL (l + e)S c (b - l)a t (k) . (23) 
n 

The halo-matter cross-power spectrum is unchanged from the /nl case, 

P hm (k, z) = (b+ 2 ^ S f-V ) Pmm{k> z) (24) 



but the halo-halo power spectrum is now 



a(k, z) 

= b f ^(k)P mm (k,z) (25) 



25 

P hh (k, z) = b}^P mn (k, z) + 4(-r JVL - fL)S 2 c (b - l) 2 P$$(fc) ■ (26) 
The second term above represents stochasticity of the halo field with respect to the dark matter field [37j . 

4. Impact of non-Gaussian initial conditions on large-scale structure: comparison 
with simulations 

To study the halo mass function and clustering, we performed collisionless A^-body simulations using the 
GADGET-2 TreePM code [Hj]. Simulations were done using periodic box size i?box = 1600 h~ x Mpc, particle 
count N p = 1024 3 , and force softening length R s = 0.05(7?box/^Vp / ' 3 )- With these parameters and the fiducial 
cosmology from ij3l the particle mass is m p = 2.92 x 10 11 h~ l Mq. Non-Gaussian initial conditions were 
implemented by generating Gaussian fields and applying the maps in Eq. Eq ([S]), or Eq. ©. The non- 
Gaussian fields were linearly evolved to the initial simulation rcdshift, z = 1 00, using the transfer functions 
from C AMB [39] . Halos were identified using the Friends of Friends algorithm [4(| with linking length L p p = 

—1/3 

0.2Ri, ox Np and halo positions are identified using the mean of the particle positions. For further details see 



In Fig. (JTJ) and Fig. ^ the ratio of the non-Gaussian mass function unc(M) to Gaussian mass function 
nc(M) is plotted. The effects of primordial non-Gaussianity are clearly visible, positive (negative) /nl, 3nl 
increase (decrease) the abundance of halos, most significantly at high mass and high redshift. The curves show 
the analytic predictions for the Edgeworth and log-Edgeworth mass functions described in £13.11 In regions where 
the effects of non-Gaussianity are small, both are in reasonable agreement with the N-body results. However, at 
the high mass end where non-Gaussian effects are largest the log-Edgeworth mass function is clearly in better 
agreement. 

In panel (b) of Fig. @ the non-Gaussian corrections for (/nl, tjvl = (|/nl) 2 ), (/nl = 0, ^nl) and (/nl, 
tml — 2(1 /nl) 2 ), are plotted together. For tml ^ (I/nl) 2 , the non-Gaussian effects are clearly larger than in 



G 
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5 



10 13 




10 14 10 15 
M (h" 1 M ) 

(a) 



a 2 



'sims, f^Uoo'T^S^f^""' 
5 

Edgeworth, f NL =±500, t nl =2 (6 fj« 

5 



log Edge., f m =±500, t nl =2 (|_f NL )= 




I I I I U I I L 



10 13 10 14 10 15 

M (h" 1 M ) 

(b) 



Figure 1: The non-Gaussian correction to the halo mass function. Plotted is ung(M, z)/ug(M, z) for (a) /nl = ±500, 
t nl = (f 500) 2 , (b) / NL = ±500, r NL = 2(| 500) 2 




10 13 10 14 10 15 10 12 10 13 10 14 10 15 

M (hr 1 M ) M (h" 1 M ) 

(a) (b) 



Figure 2: The non-Gaussian correction to the halo mass function. Plotted is (a) /nl = 0, <7nl = ±5 x 10 6 . In panel 
(b) the non-Gaussian corrections for (/nl,tjvl) = (500, (|500) 2 ), (/nl,tjvl) = (500, 2(|500) 2 ) and g N L = 5 x 10 6 are 
plotted together for comparison. 



the /NL-only case. However, the mass and redshift dependence of the curves is similar: we found that models 
with tnl > (|/nl) 2 can be made to look like models with tml = (|/nl) 2 by using a larger value of /nl- 
On the other hand, the mass dependence of (7nl effects is distinctly steeper at high masses, thus in principle 
primordial skewness and kurtosis may be distinguishable with the mass function. 

In Fig. ([3]) we show the scale-dependent bias from (?nl initial conditions. The fc-dependence of the bias is 
accurately described by the form b + gNhf3 g /ct(k) where j3 g is a constant. In panel (b) we compare the peak- 
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Figure 3: (a) An example to illustrate that halo bias in a <?nl cosmology takes the functional form form b(k) = bi + 
/3g<?NL/a(fc)- This figure corresponds to redshift 2 = 0.5 and halo mass range 1.15 < M < 1.83 x 10 14 /i -1 Mq, but we 
find the same functional form for all redshifts and halo masses, (b) The mass dependence of the <;nl bias, compared 
with the peak-background split prediction 3 f ^° " ; there is excellent agreement. 
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Figure 4: (a) An example to illustrate that halo bias in the two-field model of £12.31 with tnl = (1 + £ 2 )(§/nl) 2 is 
stochastic (ry 7^ 0). This is in contrast to the Gaussian case and the /NL-only model given in H2.ll for which halo 
bias is non-stochastic on large-scales (rij — 0). This figure corresponds to redshift z = 0.5 and halo mass range 
1.15 < M < 1.83 x 10 14 h~ x Mq, but we find the same functional form for all redshifts and halo masses, (b) An example 
of the peak-background split prediction (solid curve) with the measured stochasticity - in this case the lowest-order 
peak-background split clearly over-predicts the stochasticity. 



background split prediction (3 g = 3<7NL<9(logn)/<9/NL with the value of j3 g measured from simulations. The 
agreement is excellent. Also shown is the analytic prediction for 3<9(logn)/<9/NL using the Edgeworth mass 
function. At these masses the analytic form is not sufficiently accurate to describe the ^nl bias. 

We now consider the stochasticity induced by initial conditions with tjvl ^ (|/nl) 2 - The stochasticity 
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parameter is given by, 

Pmm(k) Pmm(k) 

where hi stands for halos in the i th mass bin. Stochasticity is generally expected on small scales from the 1-halo 
term in the halo power spectrum. But on large scales halos are expected to be non-stochastic tracers of the 
dark matter. In Fig. (Q} we show that the stochasticity vanishes on large-scales for Gaussian initial conditions 
and /nl initial conditions described by Eq. @. On the other hand, the two field initial conditions described in 
Eq. ([9]) giving tml / (f /nl) 2 do give rise to large-scale stochasticity. 

In panel (b), we show the scale dependence of the stochasticity in comparison with the lowest-order peak- 
background split prediction from H3.3I Here we find disagreement: the scale dependence of the analytic predic- 
tion agrees well with what is seen in simulations but the amplitude is too high by ~ 50%. The disagreement 
between the amplitude of the stochasticity seen in simulations and the prediction from Eq. (|26|) varies with mass 
and redshift, but is typically ~ 30% [l4j |. For Gaussian initial conditions, we also found inconsistent agreement 
between the halo model predictions for ru and the measured values. This mismatch in r,j values is in qualitative 
agreement with [4lj . 



5. Summary 

We have considered the impact of three types of primordial non-Gaussianity, described by parameters /nl, 
<7nl and tnl, on the abundance and clustering of dark matter halos. Analytic predictions for the halo abundance 
that are based on using the Edgeworth series for the non-Gaussian PDF in the Press-Schechter model agree 
well with simulations, provided the "log-Edgeworth" truncation is used (see Figs. (P), ([2]) & [13]). We found a 
simple peak-background split description of halo bias from <7nl initial conditions (Eq. (|18j)) that is in excellent 
agreement with simulations (see Fig. (J3J & [Hj]). Our simulations show that the two-field model given in ^2.31 
that gives rise to tml ^ (§/nl) 2 , does indeed generate large scale stochasticity. Unfortunately, the analytic 
description in H3.3I predicts only the fc-dependence of the stochasticity accurately, the amplitude is incorrect at 
the < 50% level (see Fig. (0} & 14]). We conclude that scale-dependent bias from /nl and <7nl initial conditions 
is well-understood analytically. On the other hand, modeling halo stochasticity appears to be more difficult and 
the analytic prescription given in ^3.31 is insufficiently accurate to interpret data without input from simulations 

0. 
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